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The 1-loop corrected decay widths of sparticles (charginos, neutralinos, gluino and sfermions) in 
the framework of the MSSM are calculated systematically using GRACE/SUSY-loop, which is 

OO ■ 

the program package for the automatic calculation of the MSSM amplitudes in the 1-loop order. 
We present the renormalization scheme used in our system and show some numerical results of 
decay widths of sfermions and gluino using the SPSla' parameter set and other SUSY parameter 
sets. 
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1. Introduction 

The supersymmetric (SUSY) theory is a good candidate for the theory beyond the standard 
model. So the experimental confirmation of the SUSY theory is one of the most important themes 
of the present and future particle experiments. Actually, discoveries of the SUSY particles (spar- 
ticles) are expected at LHC and ILC. Since we expect very accurate experimental data at ILC, 
we need theoretical prediction which match the measurement accuracy. Apparently, the tree-level 
calculation is insufficient. So we are calculating the radiative correction to possible major decay 
modes of sparticles using GRACE/SUSY-loop [j|]. In this paper we report numerical results on 
the decays channels for squarks and gluino ^ B], which are SUSY partners of quarks and gluon, 
respectively. 

2. GRACE/SUSY-loop 

Analytical evaluation of physical amplitudes characterized by many Feynman diagrams is not 
easy. It is principally for this reason that we have developed the GRACE system m\. This sys- 
tem calculates the cross sections and the decay widths and generates events automatically in the 
following way, (1) it generates all Feynman diagrams automatically, (2) it generates physical am- 
plitudes automatically, (3) it incorporates libraries (loop integral, kinematics, etc.), (4) it integrates 
the matrix element by the adaptive Monte Carlo method, (5) it generates Monte Carlo events, 
(6) it has various self-test mechanisms of the results (UV and IR cancellation, NLG invariance 
etc.). GRACE/SUSY-loop can calculate the SUSY amplitudes up to 1-loop order. For the cal- 
culation of the SUSY amplitudes at 1-loop level, there are also other programs, SloopS [§] and 
FeynArt/Calc ®. 

In GRACE/SUSY-loop, we have used the technique of the non-linear gauge (NLG) in or- 
der to test the system. Concretely, we introduce the following gauge fixing terms in the Lagrangian. 

F w ± = (d^±ieaA^±igcosd w j5Z^)W ± v±i% w ^(v + 8 H H + 8 h h ±ikG )G ± , (1) 

F Z = + ^ (v + ~e H H° + ~e h h Q )G°, (2) 
F y = d ll A». (3) 

They contain seven independent NLG-parameters, (a, j8, 8h, 5\, k, En, Eh)- We emphasize that the 
NLG interactions are included in the electroweak (ELWK) sector of the MSSM Lagrangian. While 
each Feynman diagram depends on the NLG-parameters, the sum of all diagrams should be inde- 
pendent of them. We can confirm the validity of calculation when the physical quantities do not 
change the value for different sets of numerical values of the NLG-parameters. It is the test of the 
gauge invariance. 

We use the on-mass-shell conditions as much as possible for the renormalization of the ELWK 
sector. As a result, gauge bosons, all fermions /, sfermions / and the lightest neutralino %\ have 
no mass shifts in the ELWK 1-loop order. We should note that there are some freedom in the renor- 
malization scheme of the sfermion sector. They are distinguished by different choice of residue 
conditions, decoupling conditions on the transition terms between lighter and heavier sfermions, 
and the left-handed SU(2) relations in the 1-loop order. In this paper the calculation is done with 
the scheme in which we impose the residue conditions on all sfermions except for heavier stop and 
sbottom (?2, b-i). The external line corrections for these particles become non-zero in this scheme. 
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Renormalization in the QCD sector is done in a mixed scheme. Light (1st and 2nd generation) 
quarks and gluon are treated in the DR scheme as in the convensional perturbative QCD. Massive 
particles are handled by the on-mass-shell scheme as in the ELWK sector. For the regularization 
of infrared divergences, the previous version [[[]] of the GRACE/ SUSY system used the fictitious 
mass of gluon A. We have developed a new system in which mass-singularities are regularized 
by the dimensional method. In order to refer the ultraviolet and infrared divergences we define the 
notations of Cuv = 1/e and Cm = 1/e, where the dimension of the space-time d = 4 — 2e =4 + 2e. 
In the following numerical calculation, we mainly use the SPSla' parameter set 

3. Squarks and gluino decays 

Possible decay modes of sfermions and gluino are as follows. 

q^qxh q^q'Xk (i = 1 ~ 4, k = i,2) (4) 

J^ixf, (i=l~4, * = 1,2), (5) 

g^qqj- (7=1,2) (6) 

Note that squarks cannot decay into q + g because the gluino mass is larger than all squark masses 
in the SPSla'. Gluino decays into most of quark— squark pairs, but it cannot decay into the top and 
the heavier stop ?2 because of trig <m t + mj 2 . 

Among various squark decay channels here we focus on the lighter stop t\ decays. From the 
tree level calculation, we find Brfji -»• bxt) = 86.7% and Br(t x -> t%° y ) = 13.3%. In the following 
we show detailed result on the main mode t\ — > b%± . Twelve Feynman diagrams for the 1-loop 
electro weak correction among 82 diagrams are shown in Figure 1. 




Figure 1: Twelve 1-loop electroweak Feynman diagrams for t\ — > b%l among 82 diagrams. The 82th 
diagram is counter term. 

Since GRACE is an automatic calculation system, it is important to test reliability of the results. 
The numerical values obtained by the system must pass all following tests. First, we show the inde- 
pendence of the NLG-parameters. If we change 7 gauge parameters, each evaluated value of loop 
graphs changes but the total does not change. We compare the case 1 : ( a, j8, k, £//,£/,) = 
(0,0,0,0,0,0,0) with the case 2 : (1000,2000,3000,4000,5000,6000,7000). We have obtained 
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the sum of the loop contribution, 8Ti oop and soft-photon contribution, 8T so f t as follows (units in 
GeV). 

easel : 8r hop + 8T soft = 0.15117115752797127186610833503954323 
case2 : 8Ti oop + 8T soft = 0.15117115752797127186610833480863836 

Both numbers agree up to 26 digits. So, the gauge invariance is confirmed. Next, we show cancel- 
lation tests of ultraviolet (UV) and Infrared (IR) divergence. We keep the UV divergent factor Cuv 
and a tiny fictitious mass X of photon to regularize the IR divergence in the program. In Table 1, we 
find that 8Ti oop is the same for Cuv = and Cuv = 1000. As for the cancellation test of infrared 
divergence, we see the sum of 8Ti oop and 8T so f t remain unchanged when we change X = 10~ 24 to 
10~ 27 (GeV). Finally, we can test that the sum of 8Ti oop , 8T so f t and the hard-photon contribution, 
8Thard is independent of the soft photon cut parameter k c . Since the numerical values of the sum 
8Telwk (= 8Fi oop + 8F so f t + SFhard) are the same within accuracy, the reliable value, 13.9% 
correction is obtained. 



Cuv 





1000 








X (GeV) 


io- 24 


io- 24 


io- 27 


io- 24 


k c (GeV) 


i(r 3 


10~ 3 


10~ 3 


10~ 5 


8T loop (GeV) 


-0.06256 


-0.06256 


-0.09364 


-0.06256 


8T soft (GeV) 


0.21373 


0.21373 


0.24481 


0.19301 


8T hard (GeV) 


0.04849 


0.04849 


0.04849 


0.06921 


8Telwk (GeV) 


0.19966 


0.19966 


0.19966 


0.19966 


STELWK/^tree 


13.9% 


13.9% 


13.9% 


13.9% 



Table 1: ELWK 1-loop corrections to the t\ b%+ decay width for SPSla' 



For the QCD correction (see Table 2), we use the DR scheme. Like the ultraviolet divergence 
factor Cuv, we keep the infrared divergent factor Cm in the calculation. Comparing various case of 
(Cuv,CjR,k c ), we obtain consistent and reliable value, —7.1% correction. In the GRACE system, 
we can also calculate the correction with the fictitious gluon mass and obtain the same correction 
-7.1%. 



Cuv 





1 








Cm 








1 





fec(GeV) 


IO" 3 


IO" 3 


IO" 3 


io- 4 


8F loop (GeV) 


-1.254 


-1.254 


-1.479 


-1.254 


8r soft (GeV) 


-3.752 


-3.752 


-3.527 


-4.786 


8Y hard (GeV) 


4.905 


4.905 


4.905 


5.939 


5r ecD (GeV) 


-0.100 


-0.100 


-0.100 


-0.099 


8T qcd /T tree 


-7.1% 


-7.1% 


-7.1% 


-7.1% 



Table 2: QCD 1-loop corrections to the t\ -> decay width for SPSla' 

Adding the ELWK and QCD corrections, we obtain 8Y (h b%t) l^tree = 13.9% -7.1% = 
6.8%, where T tree = 1.43GeV. Similarly we obtain the results for the other channel, 8T (?i — > t%]) 
jT tree = 7.6% +2.1% = 9.7%, where T tree = 0.22GeV. Fig. 2 shows the branching ratio at the tree 
and the 1-loop level. There is only a small difference between them. 
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Figure 2: Branching ratios of t\ decay. Outer circle 
and inner one correspond to tree level and 1-loop cor- 
rected, respectively. 
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Figure 3: nig dependence of 8 p gCT (t\ 



■bW+Xi) 



Up to now, we have considered the SPS la' parameter set. Now we take another MSSM param- 
eter set in which two-body decay channels t\ — > t%\ and t\ — \ b%\ of the lighter stop t\ is kinemati- 
cally forbidden and t\ dominantly decays into 3-body channel, t\ — > bW + %®. This is the case when 
the following mass relations hold, mj >nib+ m™ + m^o ,my <m t +m~a and m? <mh + . As an 
example, here we take = 300GeV, m^o = 195GeV, m%+ = 396GeV and obtain T tree = 0.664keV. 
Through the gluino— squark loop contributions, the corrected width depends on the gluino mass as 
well as masses of the 1st and 2nd generation squarks. In Fig. 3 the gluino mass dependence of the 
correction r QLD are shown. We find that 2TeV gluino mass shift induces about 5% shift of the 

1 tree 

width. If the gluino g is too heavy to be produced at future colliders, the precision measurements 
of the decay width of the light stop t\ will give us information of the gluino. 

Similarly, we have calculated the 1-loop ELWK and QCD correction of gluino decays (||). As 
for the QCD correction, our results agree with the calculation done by Beenakker et al. [gj], if we 
use the same input values as theirs (see Fig. 4). When we take the SPS la' parameter set, we obtain 



* r(ref) 1Joop 

r(GRACE) Iree 

r(GRACE) 1Joop 



550 600 

m„iui„„ (GeV) 



Figure 4: Gluino mass dependence of QCD correction to tt\ . Dotted line and solid line respectively cor- 
responds to r tree and ri_i 00 p calculated by GRACE. Rectangle and triangle corresponds to Ttree and Ti_i 00 p 
shown in Fig. 5 of the reference [BR. 



8T/T tree (g^bbi) = -18.0%(QCD) + 2.5%(ELWK) 
8T/T tree {g-> tii) = -13.0%(QCD) + 1.2%(ELWK). 



(7) 
(8) 
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4. Summary 

Since the sparticles are expected to have masses of the order of the ELWK scale, we cannot ne- 
glect the ELWK corrections as well as the QCD corrections in the precise theoretical prediction of 
their production cross sections and the decay rates. Using GRACE/SUSY-loop, we have system- 
atically calculated ELWK and QCD corrections to the sfermion and gluino decays, and confirmed 
that both corrections are equally important. We have found that QCD corrections for t\ — > is 
the same in the DR regularization and in the regularization with the fictitious gluon mass. We have 
already calculated the chargino [|T]] and the neutralino decay channels and have got reliable results. 
Extension of the adaptive range to radiative corrections for multi-body channels is now planned. 

This work is partially supported by Grant-in-Aid for Scientific Research(B) (20340063) and 
Grant-in-Aid for Scientific Research on Innovative Areas (21105513). 
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